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ABSTRACT 

An ideal pulse tube cryocooler using an ideal gas can operate at any temperature. 
This is not true for real gasses. The enthalpy flow resulting from the real gas effects of 
He, 4 He, and their mixtures in ideal pulse tube cryocoolers puts limits on the operating 
temperature of pulse tube cryocoolers. The discussion of these effects follows a previous 
description of the real gas effects in ideal pulse tube cryocoolers and makes use of models 
of the thermophysical properties of 3 He and 4 He. Published data is used to extend the 
analysis to mixtures of 3 He and 4 He. 

The analysis was done for pressures below 2 MPa and temperatures below 2.5 K. 
Both gasses and their mixtures show low temperature limits for pulse tube cryocoolers. 
These limits are in the 0. 5-2.2 K range and depend on pressure and mixture. In some 
circumstances, even lower temperatures may be possible. Pulse tube cryocoolers using the 
two-fluid properties of dilute 3 He in superfluid 4 He appear to have no limit. 
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INTRODUCTION 

An ideal pulse tube cryocooler using an ideal gas has a cooling power given by [1] 

Q c /W = T c /T h (1) 

Such a cooler does not have a minimum operating temperature. It is capable of cooling to 
any temperature. When real gasses are used, the performance of pulse tube cryocoolers is 
affected at all temperatures [2-6]. The performance is significantly reduced for 
temperatures near and below the critical temperature of the working fluid. This was shown 
analytically for 4 He [2-3] and 3 He [5]. Since the critical point of 3 He occurs at a lower 
temperature than for 4 He, pulse tube coolers using 3 He reach lower temperatures than those 
using 4 He. This has been demonstrated in practice [3-4], The existence of low temperature 
limits of pulse tube cryocoolers were introduced in a recent study [7], 



This paper reports the results of an analysis of the ultimate operating temperatures for 
an ideal pulse tube cryocooler for pressures below 2 MPa using 3 He, 4 He, and their 
mixtures as the working fluid in ideal pulse tube cryocoolers. 


APPROACH 


The approach has been previously described for ideal and non-ideal pulse tubes [8-9], 
for 4 He real gas effects in ideal pulse tube cryocoolers [1], and for 3 He real gas effects in 
ideal pulse tube cryocoolers [5]. 

We make the following assumptions. Pressure gradients that result from viscous 
flows in the orifice or inertance tube are proportional to and in phase with the mass flow. 
Similarly, heat transfer is proportional to and in phase with temperature differences. All 
oscillating quantities are sinusoidal at the fundamental frequency. These are the 
assumptions of small amplitudes and that the energy carried by higher harmonics is small 
compared to the total energy carried by the system and can be ignored. 

The approach assumes that the system (FIGURE 1) is ideal except for the gas 
properties, which are the real gas properties. The gas properties for pure 3 He and 4 He were 
taken from thermophysical models [10], In ideal, lossless components, there is no entropy 
generation. Thus, the effects of viscous flow losses (except in the orifice or inertance 
tube), secondary flows, and inadequate regenerator heat capacity are ignored. Heat flows 
(thermal conduction) are treated as an external load. Only steady state operation is 
considered. The system control volume is a closed system (no mass flows in or out). 

In steady state operation the mean enthalpy and entropy flows are given by 


X = 



mxdt 


where x is the period of oscillation, x is h or 5 and X is H or S, respectively. 


( 2 ) 


IDEAL GAS PROPERTIES 

If the working fluid is assumed to be an ideal gas, then 



ds = cp dT/T- R dP/P 

( 3 ) 

and 

dh = cp dT 

( 4 ) 
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FIGURE 1. The principal components of an Orifice Pulse Tube Cryocooler showing the location and 
nominal direction of the heat and work for a control volume enclosing the entire cryocooler. 




where cp = 5/2 R for helium. 

In the pulse tube, an adiabatic region, ds = 0, thus equations (3-4) become 



dT/dP = V/c P 

(5) 

and 

dh = V dP 

(6) 


where the coefficient of dP in equation (6) is evaluated at the temperature of the cold heat 
exchanger. Elsewhere in the pulse tube, VP increases with temperature. However, because 
no work is done and the pulse tube is adiabatic, there is no change in the enthalpy flow 
within the pulse tube. To maintain constant enthalpy flow, the mass flow decreases with 
increasing temperature. This accumulation effect was discussed earlier for ideal gas flows 
[9]. In adiabatic sections the temperature and pressure oscillations are in phase. 

In the regenerator, there are no temperature oscillations, dT = 0, thus equations (3-4) 
become 



ds = (V/T) dP 

(7) 

and 

dh = 0. 

(8) 


Thus there is no enthalpy flow in the regenerator. 


REAL GAS PROPERTIES 

When real gas effects are included, equations (3-4) become [11] 

T ds = cp dT - T ft V dP (9) 

and dh = cp dT + [1 - T$]VdP (10) 

where (3 = \/V dV/dT\ P is the volume expansivity. Also, c P is no longer 5/2 R, but is a 
function of T and P. For an ideal gas, 7|3 = 1. 

Pulse Tube 


Within the pulse tube, the flow is adiabatic. Thus, ds = 0 and equation (9-10) become 



dT/dP =7(3 V/cp 

(11) 

and 

dh = VdP . 

(12) 


where, as in equation (6), the coefficient of dP in equation (12) is evaluated at the 
temperature of the cold heat exchanger. 

Regenerator 

In regenerators, heat transfer to the regenerator material suppresses the gas 
temperature amplitude. Thus, dT= 0. For a real gas, equations (9-10) become: 



ds = - V (3 dP 


(13) 


and 


dh = [\-T$\ VdP. 


(14) 


DISCUSSION 

The heat absorbed at the cold heat exchanger is given by the difference at the cold 
heat exchanger between the enthalpy flows in the regenerator and in the pulse tube. For an 
ideal gas, this is the difference between equations (5) and (7): 

A dh=VdP. (15) 

where equations (6) and (15) have been combined. 

For a real gas, the heat absorbed is the difference between equations (12) and (13): 

A dh = T$VdP (16) 


Pure 3 He and 4 He Fluids 

Equation (15) places no limits on the temperature at which an ideal pulse tube cooler 
can operate using an ideal gas. Equation (16) does place temperature limits on an ideal 
pulse tube cooler using a real gas. The cooling reaches a limit when (3 = 0 and cooling 
ceases. The locus of (3 = 0 are shown in FIGURE 2. The 3 Fle curve and the 4 He curve near 
2.2 K were calculated from reference [10], The lower 4 He curve was calculated using 
reference [12], For 4 He, (3 < 0 between the two lines marked 4 He and (3 > 0 elsewhere. For 
3 He, (3 < 0 for temperatures below the line and (3 > 0 for higher temperatures. 

In the regions where (3 > 0, pulse tube coolers will operate in the usual fashion with 
the mass flow and pressure oscillations close to being in phase [13]. For coolers starting 
near room temperature, a cooler could reach - 2.2 K using 4 He and- 1 K using ’He. 



Figure 2. The locus in pressure-temperature space of where the cooling power of an ideal pulse tube 
cryocooler is zero for ’He and 4 He working fluids. Also shown is the lambda-line of 4 He. 


Below ~ 1 K, the 4 He (3 is again > 0. A cooler operating in this region would need to reject 
heat to a pulse tube cooler using 'He at ~ IK. Unfortunately, the 4 He (3 is very small 
below 1 K (see FIGURE 3), making it difficult to realize a practical cooler in this region. 

In the regions where (3 < 0, a cooler using the usual phase shift between mass flow 
and pressure would produce negative cooling (heating). Perhaps one could reverse the 
roles of the pulse tube heat exchangers to make a cooler; or perhaps one could change the 
phase between mass flow and pressure by an additional 180°. In practice, neither approach 
is likely to work. The 3 He (3 is quite small below 1 K (see FIGURE 3) and 4 He has a high 
thermal conductivity below the lambda-transition that will thermally short the pulse tube. 

Mixed 3 He and 4 He Fluids 

Mixtures of 3 He and 4 He offer the possibility of adjusting the gas properties to modify 
performance. While the thermodynamic state of a pure fluid can be specified by two 
parameters such as T and P, a two-component mixture requires additional parameters, xj, 
the molar fraction of 3 He in the mixture and the chemical potential, p, of each isotope. In 
steady state operation, a pulse tube cooler will have concentration gradients caused by the 
two isotopes having different viscosities and by chemical potential effects. The ideal pulse 
tube cooler model ignores viscous effects except for the pressure drop in the orifice. We 
will assume that above the lambda-line that the concentration gradients do not significantly 
contribute to the ultimate temperature and the ultimate temperature is determined only by 
the concentration at the cold heat exchanger. I.e., we assume that the loci of dV/dT\p ^ = 0 
and dV/dT\px 3 = 0 are the same. 

There is only limited data on the volume expansivity of helium mixtures in the 1 -2 K 
and < 2 MPa temperature and pressure range. In FIGURE 4, (3 = 0 data are shown as a 
function of mixture ratio superimposed on the 3 He/ 4 He phase diagram. All of the data are 
for saturated vapor pressure rather than constant pressure. The mixture data is for constant 

X3. 

The pure fluid data points for the upper branch of (3 = 0 (JC3 = 0 and X3 = 1) were 
calculated [10]. The other two data points are interpolated from a table [14], In addition to 
errors in the published data, the accuracy of the interpolation is about ± 40 mK because of 




Figure 3. Volume expansivity of 3 He (left) and 4 He (right). 



Figure 4. The (3 = 0 data for 3 He/ 4 He mixtures at saturated vapor pressure superimposed on the 3 He/ 4 He phase 
diagram. The data are discussed in the text. 

the large values of the second and higher order derivatives near (3 = 0. These four data 
points are well fit by a linear function of x 3 , thus, further justifying the assumption of non- 
interacting components. All the upper branch of (3 = 0 lies in the He-I region. 

The pure fluid data point for the lower branch of (3 = 0 (x 3 = 0) was calculated [12], 
The other data point is an upper limit estimated from data [15]. The actual (3 > 0 region is 
likely smaller than indicated. 

3 He Quasi Particles in Superfluid 4 He 

The proceeding discussion assumes that the mixture moves as a bulk fluid. In the 
superfluid region, He-11, the fluid can act as two interpenetrating fluids, which can move 
independently. One is the pure 4 He superfluid component. The other is the normal fluid 
component consisting of some 4 He and all of the 'He. Near the X-line, most of the fluid is 
in the normal state. Below ~ 1 K, the normal fluid consists of only 3 He quasi particles. If 
the piston in the compressor has a superleak then the He can act as a low-density gas 
moving in a fixed 4 He background. A full treatment of the thermodynamics in this regime 
is beyond the scope of this paper. Reviews of the thermodynamics of dilute 3 He in 4 He 
solutions have been given elsewhere [16-18]. 

A pulse tube using dilute 'He in a fixed 4 He background was developed at Los 
Alamos. It was filled with x 3 = 0. 17 and operated between 1 K and 0.6 K [19]. The system 
pressure is fixed by the conditions in the compressor/aftercooler. The compressor does not 
cause pressure oscillations; rather it causes the 3 He concentration, X 3 , and the osmotic 
pressure, n 3 , to oscillate. Pressure gradients in the regenerator, heat exchangers, and 
orifices are replaced by gradients in n 3 . In the pulse tube, the constraint that VP = 0 is 
replaced by Vq 4 = 0, where m is the chemical potential of the 4 He. Several loci of Vp 4 = 0 
(from data in reference [16]) are shown in FIGURE 5 from which one sees that in the pulse 
tube X 3 increases as the temperature decreases. The cold end temperature is limited by 
phase separation. To reach the lowest temperatures limits x 3 < 0.064. Such a pulse tube 
with a warm end above 0.6 K will have x 3 < 0.01 at the warm end and in the compressor. 
Such low concentrations limits the mass flow and the cooling achievable in a practical 



Figure 5. Selected loci of V[x 4 = 0 for dilute 'He/ 4 He mixtures. The dashed line is, approximately, the 
operating region of the Los Alamos pulse tube. 

sized cooler. This type of cooler has not seen further development nor has it seen any 
practical application 

The present thermodynamic model of dilute 3 He in 4 He [16-18] shows no evidence of 
a low temperature limit of a pulse tube cryocooler. 


SUMMARY 

The theoretical limit of a non-superfluid pulse tube cryocooler is when the volume 
expansivity is zero: (3 = 0. Below ~ 1 K, 3 He concentration driven pulse tube cryocoolers 
are possible with no known ultimate limiting temperature. 

In practice, a conventional pulse tube cryocooler using 4 He is limited to temperatures 
above ~ 2.2 K. Using 3 He, the limit is above ~ 1 K. The limit can be adjusted in between 
these two points by using a mixture of 4 He and ’He. When mixtures are used, the mixture 
ratio will not be constant throughout the cooler. The mixture ratio will vary with position 
because the chemical potential of the mixture depends on temperature and pressure. The 
mixture at the cold heat exchanger will not be the same as the initial mixture; rather, the 
mixture stratifies. Thus, the ultimate achievable temperature of a cooler using a mixture 
will depend on the mixture ratio at the cold heat exchanger when the cooler is operating in 
steady state. 
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NOMENCLATURE 


Roman 

H Enthalpy 

P Pressure 

Q Heat 

R Gas constant (per mole) 

S Entropy 

T Temperature 

V Specific volume 

W Work 


c Specific heat 

h Specific enthalpy 

m Mass 

i Specific entropy 

x Molar fraction of ’He in 

mixture 


Greek 

|3 Thermal expansivity 

M Chemical potential 

n Osmotic pressure 

Subscripts 

3 3 He component 

c Cold 

h Hot 

P Isobaric 
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